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We calculated the differential and total rates for the semileptonic decays of K1 and mesons 
taking into account a linear dependence of the Ki^ form factors. As a case in point, we included 
these rates into the calculation of the atmospheric neutrino flux at energies 1 to 100 TeV. The 
calculated neutrino spectra are between the earlier predictions while the neutrino flavor ratio, R, is 
somewhat affected by the iffs form factors. The R proves to be very sensitive to the contribution 
, of neutrinos from decay of charmed particles, providing an additional method to test the charm 

■ production models in future experiments with large-volume neutrino telescopes. 

^ ' PACS numbers; 13.20.Eb, 13.85.Tp, 96.40.Tv 

'. I. INTRODUCTION 

cn ; 

' Atmospheric neutrinos (AN) come from the decays of unstable particles generated in the collisions of primary 
and secondary cosmic rays with air nuclei. Up to very high energies, the dominant contribution is due to decays of 
charged pions, charged and long-lived neutral kaons. The AN from these sources have come to be knovifn as "tt, K" 
^ or conventional neutrinos. As energy increases, semileptonic decays of charmed hadrons (mainly , D'^, D^, Df^ 
^—1 . mesons and A+ hyperons) should become important. These AN are usually called prompt neutrinos (PN). The 

■ borderline energy between the domination regions of conventional and prompt neutrinos is a long-standing question. 
C^l Taking into account the available data on cosmic-ray muons [Q, it is safe to say that the PN fraction in the vertical 

(horizontal) muon neutrino flux is negligible at energies below E"^ ~ 1 TeV (~ 10 TeV) ; for the electron neutrino flux, 
2^ , the corresponding borders are roughly an order of magnitude less. But it is not excluded that may be increased 
^ ■ by 10 to 100 times. 

[ The AN flux with energies above 1 TeV represents an unavoidable background for many of the astrophysical 
JpH, experiments with the full-size underwater/ice neutrino telescope^. At the same time, the AN flux is a natural tool 
O ^ for studying neutrino interactions with matter and, along with the atmospheric muon flux, it provides a way of testing 
the inputs of nuclear cascade models that is parameters of the primary cosmic-ray flux and cross sections of hadron- 
nucleus interactions at energies beyond the reach of accelerator experiments. In particular, the AN flux measurements 
J> , have much potential for yielding information about the mechanism of charm production. 

In both cases - to correct for the AN background and to use the AN flux as an additional tool of particle physics 
- it is necessary to know with confidence the conventional AN flux. The relevant calculations have been performed 
by many authors (see [^-|8| and references therein). The discrepancy between the different predictions for the tt,K 
neutrino flux above 1 TeV ranges up to 25-35 % for -I- 77^ and to 60-70 % for i^e -I- (depending upon the zenith 
angle). Unfortunately, these numbers are not representative of the upper limits for the calculation uncertainty. For 
the most part this is due to the incompleteness in the current knowledge of the primary spectrum and composition 
as well as mechanisms of tt and, to a greater extent, K meson production at high energies (see |^ for a detailed 
discussion). A sizable part of the discrepancy is caused by different approximations and simplifications employed in 
the calculations. 

Predictions for the PN contribution vary by a few orders of magnitude (see and references therein). Here, the 
basic challenge is of course in the mechanism of charm hadroproduction. An additional (though not so drastic) source 
of uncertainty has to do with the differential rates of inclusive semileptonic decays of charmed hadrons. The theory 
of charm production and decay is still far from completion and the corresponding accelerator data are rather meagre. 
However, when correlating the predictions of different models with each other and with experiment, the neglecting 
some comparatively small effects may introduce additional systematic errors and should be avoided. 
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^Among them the detection of neutrinos from the (quasi)diffuse neutrino backgrounds, like pregalactic neutrinos, neutrinos 
from the bright phase of galaxy evolution and from active galactic nuclei. 
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Some relative characteristics of the AN flux, like the zenith-angle distribution and the flavor ratio, 



R 



prove to be less sensitive to the uncertainties of the conventional AN flux predictions and, on the other hand, they are 
very dependent of the PN contribution. The latter is owing to the salient and model-independent features of the PN 
flux. First, it is almost isotropic within a wide energy rangeQ and second, the v/V ratio and the flavour ratio, R, are 
both about 1. These features provide a way to discriminate the PN contribution through the analysis of the angular 
distribution and the relationship between the muon and electron neutrino-induced events in a neutrino telescope. 

In this paper we calculate the conventional AN energy spectra with taking into account one additional effect which 
was ignored in all previous AN flux calculations, - the dependence of the and K^^ decay form factors on (where 
q is the 4-momentum transfer to the lepton pair). Inclusion of the g^-dependent form factors causes certain changes 
in the differential Ki^ decay rates and therefore it should affect the overall tt, K neutrino flux. Magnitude of the effect 
is wittingly small, but it is energy-dependent and opposite in sign for muon and electron neutrinos; consequently it 
must change predictions for the i? as a function of energy. The Ki^, decay contribution is "by definition" negligible at 
energies where the prompt neutrinos become important. But, let us recall, the borderline energy remains unknown 
for the present and hence the effect under consideration might be interesting up to the PeV energy range. 

In the next Section we present some necessary formulas from the theory of neutrino production in the atmosphere. 
Section [II deals with the differential and total K^^ decay rates. In Section |^ we briefly describe our nuclear cascade 
model. Section ^ is devoted to the discussion of our results for the AN spectra, while the conclusions are in Section VI . 



II. NEUTRINO PRODUCTION IN THE ATMOSPHERE 



Let P labels a particle which can produce a lepton pair Ivg or tPi {i = e,fj,) on decay and <I>p {Ep^z,d) be the 
differential energy spectrum of these particles as a function of energy Ep, atmospheric depth z and zenith angle i?. Let 
dTp^^/dEi, be the differential with respect to neutrino energy E^, rate for the /c-body decay mode Pik, as a function 
of Ep and E^, in the laboratory frame of reference (the symbol ly stands for ve or Vi). Then the differential energy 
spectrum of neutrinos at the depth z and zenith angle d is given by 



P k 



dP^^^ {Ep,E,) 



dE^ 



$P {Ep,z',d) 



dz'dEp 



(2.1) 



where p{z^'d) is the air density on the corresponding altitude in the atmosphere in terms of variables z and -d: 
summation is over all particles P and fc-body decay modes. 



the 



The main decay modes answerable for the conventional neutrino production are ^, 



± 



e3' ^i2^ ^'■^l2^ 



± 



K^27 ^^"^ ^es- Prompt 



neutrinos are produced through the multiple modes of semileptonic decays of charmed hadrons but, considering sizable 
gaps in the experimental data |10| and certain difficulties in the theoretical description of charm decay, it is instructive 
(and generally accepted) to use the inclusive ap proa ch. 

The lower limit of integration over Ep in eq. ( |2.l|) is defined by kinematics {c — 1). At ^ 1 GeV, 



Ef 



{ml/MljE^ 



where 



= mi- — m,, 



ML. = 



and, for the inclusive decays. 



M - = m 



^Namely, at 10 TeV < < 3 ■ 10^ TeV, the maximum anisotropy is about 3-4 % 
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with s™'" the invariant mass square minimum. 

The total rate of the Pik decay in the lab. frame is defined by 



B{Pik)mp 

dE^ = , 

TpEp 



where mp and rp are the mass and life time of the particle P, respectively, and B (Pik) is the Pik decay branching 
ratio. Let us introduce the "spectral function" 



r {Pek 



dr^,^ iEp,E,) 



dE,, 



It can be shown that, in the ultrarelativistic limit, the Fp^^ is a function of the only dimensionless variable x = E^/Ep 
The spectral function for any two-body decay is merely constant. In particular. 



1 „u.,(V,.-) 1 



1- ml/ml' ^'^^ 

The spectral functions for the three-particle decay of a polarized muon in the ultrarelativistic limit are of the form p^ ] 

Kf'^i^) ^ 2(1 - [l + 2x± V^{1 - 4x)] , 

^Mei^^'H^;) = -^(1 - a;) [5 + 5a; - 4a;2 ± -p^(l + x - %x^)\ , 
o 

where "P^ is the muon polarization dependent on the muon and parent meson momenta. These dependencies are 
different for the different meson decay modes. The /x-decay contribution into the + flux is very small at high 
energies but in contrast, it dominates in the Ve + T'e flux up to about 100 GeV for vertical and to several hundreds 
of GeVs for horizontal directions. The muon polarization is therefore an essential factor affecting the neutrino flavor 
ratio and the neutrino to antineutrino ratio However, ai Ey > 1 TeV one can greatly simplify matter treating 

the Vfj, as an effective constant, (P^). In our calculations we adopt (Vf^) — 0.33. Besides, as is customary in all AN 
flux calculations, we take no account of a small change of the shape of neutrino distributions which result from the 
radiative mode fj, eVf.v^'y (with branching ratio of (1.4 ± 0.4) %) but simply increase the B (fies) to 100 %. 

The spectral functions for the three-particle kaon decays calculated without considering the if ^3 form factors can be 
found in |^ and The inclusion of the q'^-dcpendent form factors is the subject of the next Section. One more point 
need to be made here. As in the case with /j,-decay, below we neglect the radiative mode — > 7reFe7 (branching 
ratio is (1.3 ± 0.8) %) but increase the B {kI^^} to 40% (cf. ||l^). This approximation yields a completely negligible 
change in the + Ve flux. 

In this paper, we will not enlarge on the calculation of the spectral functions for the inclusive semileptonic decays 
of charmed hadrons. Our simple phenomenological approach to the problem has been described in [Q. 

III. Ki3 DECAYS 

In the standard theory of weak interactions, the Ki^ decay matrix element can be written in the form 

^ sin 0c [f+{q^){pK+P.r + f-{<l^)(pK-p.r]h,{l^l^W (3.1) 

Here Gp and Oc are the Fermi constant and Cabibbo angle, pk,tt are the 4-momenta of the kaon and pion, and f±{q^) 
are dimensionless form factors which are real functions of = {pk —Ptt)'^, the square of the 4- momentum carried by 
leptons. Experimental investigations of decays suggest that the form factors f±{q^) are smooth functions of q^ 
which are normally written in the linear approximation as 

-.2^ _ /n^ / 1 , \ . 1 



/±(9')-/±(0) 1 + A±^ . (3.2) 
In the limit of unbroken SU{i) symmetry, 

/+(0) = 1 for ifOg and /+(0) = ^ for if±, (3.3) 
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while /-(O) reduces to zero. As a consequence, the parameter ^ = /_(0)//+(0) should be small for K^s decays (the 
Ademollo-Gatto theorem) In the if^a case, the Ademollo-Gatto theorem is not valid and so, it is not forbidden 
that ^ 1. The absolute normalization of the Ke^ decay rates is n ot warr anted for our purposes, as we use the 
experimental values for B{Ke3) and tk- This being so, we adopt eqs. (3.2 3.3) from here on, considering X± and ^ as 
input parame ters . 

From eqs. (3^) and (|3.2| ), using standard techniques [Q, we find the differential (with respect to neutrino energy) 
and total Ke^ decay rates in the lab. frame of reference. In a general way, the Kes spectral function may be written 
as 111 



with the normalization constant 



1 3 



Z= J2 C„[^n+l(0)- J„(0)], 



Here 



J„(x) 



dyy"^ y/{y- y+)iy - y-), 



y± 



v+ 



'if 



(these integrals are expressible in terms of elementary functions) and C„ are the constants proportional to /+(0) 
and dependent on the masses, m/f, m^, rrii and the parameters X± and ^. Let us define C„ ~ /^(0)c„. Then the 
coefficients c„ are 

C-4 = -&ri{r - n)^v^X'^, 



c_3 = 8re{r — ri) uvX — 2(r — rf){4(r — ri)[r — (1 — 4w)] 
-r43r(l + 3r) - r^(9 + lOr - ri)]v^}X'^, 



c_2 = -3rf {r - rt)u^ + 4{2(r - rt)[2{r - re) + 3ri{u + v)] 
~-ri[r + r£ + (r - r^)(3 + 4r - ri)]uv}X 

+8{{r - ri)[{r - + 3r + n) - (r + r^)] -f \ri{r(2 + 3r) - ri[\ + 'ir)\v}X^ 
-2ri,{ir - r<,)[3(2 + 8r + 3r^) - (2r + n)] + (r + )(3 + 3r + ri)}v^X'^, 



c_i = -3{4r - re[4{l - 2u) + (1 + r - ri)u^]} - 4{8r(l -f r) - 4rf (4 + r + \ 

+ r + r^)]uu — &ri{\ + 2r)(?i + f)}A 
-8{3r(l + 3r + r^) - re,[b + 16r + + r^(l + r + r^) 
-h2(5 -I- 15r + + 3r<.)t;}A2 + 2rf{3(l + r)(l + 9r + r^) 
-hr43 + 4r + + r^(4 + r + re)]}v^X^, 



Co = 3{4(1 + r) - r£[4 + (u - 8)u]} + 4{4(1 + 4r + r^) 

-r44(5 - r + 2re) - 6{2 + r + re){u + v) + (4 + 4r + ri)uv]}X 

+8{{1 + r)(l + 8r + r^- re[3 - ir - + re{3 -r + 3re)]}X^ 

+2re{16[3 + 6r + r^ + ri{3 + r + r^)] - [9(1 + 3r + r^) + (7 + 4r + ri)]v}vX'^, 

ci = -12 - 8{4(1 + r) - re[2 - 3(u + t;) - H}A - 2{12(1 + 3r + r^) 
-r44(7 + 5re) - 32(3 + 2r + 2ri)v + (9 -f 9r -I- 5re)v^]}X^, 

C2 16A + 2{12(1 -f- r) ~ ri\\ - (16 - 3^;)^;]}A^ 



C3 = -8A' 
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where we used the following notation: 



and 



To estimate the numerical values for the C„, we use the parameters A+ and ^ evaluated by the Particle Data Group |l^] 
(see table |). 

TABLE I. Ke3 form factor parameters (lOl. 













A+ 


0.0300 ± 0.0016 


0.034 ± 0.005 


0.0286 ± 0.0022 


0.033 ± 0.008 




-0.11 ±0.09 




-0.35 ±0.15 





Since the term of the matrix element (3.1) proportional to /„ can be neglected for Kes and most if ^3 data are 
adequately described with a constant /- in the subsequent discussion we assume A_ — 0. The numerical values 
of the coefficients C„ are given in tables |ll and III . For comparison we also included in these tables the C„ calculated 
with A-|_ — 0. 

In the case, electron mass can be neglected as compared to pion mass. In this approximation, one can easily 
found 



n=0 
C-3 



il~x) 



rC-2 - C- 
1 - X 

n+l / 1 ~ 2; 



± (C_2 - rC_i) ln(l -x) + C_i(l - x) 

„n+2 



2r 



n + 2 n+l J ' {n + l){n + 2) 
C_2 (In r ± 1) ± rC_ 1 (In r - 1 



(3.4) 



where 



Z = [(1 ± r)C-2 - C_3 - rC-i] \ni 



2(l-r)C_2-^^(C_3±rC_i) 



2r 



r(l - r"+i) 



1 - r 



ri+3 



.(n±l)(n±2) (n±2)(n±3)_ 
and r = rn^/mj^. The total K^s decay rate in the rest frame has the form (cf. [^) 



(3.5) 



Gl sin^ e. 



fl{0) (ao±aiA+ ±a2A^) 



7687r3 

with constant a^. Substituting the numerical values of the parameters yields 



(3.6) 



ao 
ao 



0.576, ai 
0.560, ai 



2.140, 
1.947, 



0-2 

a2 



1.580 for K^^, 
1.345 for if°3. 



Let us pass over the explicit form of the F^' ^{x) and T*{Kfj_3) which is much more complicated in comparison with 
eqs. (|3.4|-|3.6[), and proceed to some numerical results. Figures |l| and || show respectively the Kfs spectral functions 
and differential decay rates (versus a:) calculated with /+ = /+(0) and with the g^-dependent form factors. As is seen 
from fig. |l|, the effect for the spectral functions is different in magnitude and, more importantly, opposite for electron 
and muon neutrinos, even though the absolute differential rates (fig. H) grow if the form factors are accounted for. 
As a result, the K^^ contribution to the AN flux shall slightly be increased for electron neutrinos and decreased for 
muon neutrinos, compared to the case of constant form factors. Clearly, in the range where the conventional neutrinos 
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TABLE II. Coefficients C„ for Kl^ decays (A_ = 0, ^ = -0.11). 





- 




Kl - 




A+ = 


A+ = 0.030 


A+ = 


X+ = 0.034 


C-4 











-1.99 ■ 10-** 


C_3 





-5.90 • 10"'' 





3.78 • lO-'^ 


C-2 





1.03 • 10"^ 


-2.29 • 10-* 


4.33 • 10-* 


C-l 


-3.93 • 10-2 


-6.45 • 10-2 


-5.93- 10-2 


-7.94 • 10-2 


Co 


5.39 ■ 10"^ 


7.29 ■ 10"^ 


5.60- 10-^ 


7.80 • 10-^ 


Ci 


-5.00 ■ 10^1 


-8.19- 10^^ 


-5.00- 10^^ 


-8.82 • 10"^ 


C2 





1.66- 10"^ 





1.96 • 10"^ 


C3 





-1.21 • 10-2 





-1.56- 10-2 



TABLE III. Coefficients C„ for K^^ decays (A_ = 0, C = -0.35). 



Cn 






- 




x+ = o 


A+ = 0.0286 


\+ =0 


A+ = 0.0330 


C„4 











-6.77 ■ 10-" 


C-3 





-2.55 • 10-*^ 





1.95 ■ 10-® 


C-2 





4.66 • 10-* 


-1.51 ■ 10"* 


1.17- 10-* 


C-l 


-1.87 • 10-2 


-3.06 • 10-2 


-3.28 • 10-2 


-4.31 ■ 10-2 


Co 


2.69 ■ 10-^ 


3.62 • 10-1 


2.83 ■ 10-1 


3.96- 10-1 


Ci 


-2.50 • 10-1 


-4.10 • 10-1 


-2.50-10-1 


-4.46 • 10-1 


C2 





8.34 • 10-2 





1.01 • 10-1 


C3 





-6.10 • 10-3 





-8.12 • 10-3 



TABLE IV. Ki3 decay rates. 



Decay 


Calculated with 


Calculated with 


Experimental 


mode 


/+ = /+(0) (lO'' s-i) 


/+ = /+(g2) (10® s-i) 


best fit (10*^ s-i) 




6.76 


7.49 


7.49 ± 0.06 




4.38 


5.25 


5.25 ±0.05 


/± 

^e3 


3.34 


3.70 


3.89 ±0.05 


K% 


2.06 


2.50 


2.57 ±0.06 



dominate, the magnitude of the effect in the overall AN flux is determined by the kaon production cross sections that 
is by the "^Vtt ratio" and, to a lesser extend, by the cross sections for kaon regeneration. 

In table IV we give the deca y ra tes, T*{Ki^), calculated using the above formulas with /+ = /+(0) and with the 
linear dependence of /+ on (3.2), together with the best fits of experimental data obtained by the Particle Data 
Group 0. 

As table ^ suggests, the inclusion of the g^-dependent form factors causes the increase of about 11% in the K,.^ 
decay rates and of about 20 % in the if^a rates. The improved rates correlate well with the experimental data except 
the Kf^ case. But the latter is not an essential flaw, having regard to the variance of the world data on r*(i^^) which 
is far more than the error of the PDG best flt. The additional corrections due to the SU{'i) symmetry breaking should 
be less than 3-4% and, what is more essential for our study, they cannot change the spectral functions, F'^^^{x), that 
is the shape of the neutrino distributions from if^s decays. 



IV. NUCLEAR CASCADE MODEL 

Our nuclear cascade calculations are based on the model by Vail et al. |]l6| . The results obtained within this model 
agree well with all available experimental data on hadron spectra (including the single proton, neutron and pion 
fluxes) for various atmospheric depths at energies from 1 TeV up to about 600 TeV. The muon spectrum calculated 
with this model is in reasonable agreement with the current sea-level and underground data It is also in very 
good agreement (within 5%) with the recent Monte Carlo calculation by Agrawal et al. 

The model takes into account the processes of regeneration and recharging of nucleons and pions as well as produc- 
tion of kaons, nucleons and charmed hadrons (-D^, I?", A+) in pion-nucleus collisions, muon energy loss, etc. Let 
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us briefly enumerate here the basic assumptions and discuss some pluses and minuses of the model. Further details 
and numerical values of the input parameters can be found in 1 16 for tt and K meson production and in for charm 
production. 

(i) Our (analytical) model of the all-particle primary spectra and chemical composition is taken from Nikolsky et 
al. p7| ]. This model adequately describes the world data on primary cosmic rays from about 1 TeV/nuclcon up to the 
range well beyond the "knee" region (see a discussion in Bugaev et al. The nuclear component of the primary 
flux is treated on the basis of the superposition model. 

(ii) We assume a logarithmic growth with energy of the total inelastic cross sections cr^'Ji' for interactions of a 
hadron h with a nuclear target A. Such dependence arises from a model for elastic amplitude of hadron-hadron 
collisions, based on the conception of double pomeron with supercritical intercept. For simplicity sake we also use 
another consequence of this model, the asymptotic equality of the inelastic cross sections for any projectile hadron. 
Namely, we assume 

(^ri\Eh) = ahA + (TA HEh/Ei), for h^N,T:, K, D, A„ 

{A is the "air nucleus") at Eh > Ei = 1 TeV. The calculated cross sections for h — N,tt,K are in reasonable 
agreement with available accelerator and cosmic ray data. There are no the data on the inelastic cross sections for 
charmed hadrons but we notice that the prompt lepton flux is scarcely affected by the specific values of these cross 
sections up to about 10** TeV of lepton energy (due to very short lifetime of charmed particles). Thus even a rough 
estimation of cr]^^' and cr™*^^ will suffice for our purposes. 

(Hi) We assume Feynman scaling in the fragmentation range for the inclusive processes hiA — > hfX, with hi — 
p, 71, TT^ and hf — p,n, tt^, , , . Thus the invariant inclusive cross sections Ed'^a^^/d^p are energy independent 
at large Feynman x. The truth of this assumption is an outstanding question. 

(iw) The kaon regeneration (i.e. the processes K^A — > K^X, K^A K'^X, and so on) is disregarded in our 
calculations. We also neglect nucleon, pion and charm production in kaon-nucleus collisions as well as pion production 
in kaon decays, which makes it possible to study the "ttTV" cascade without reference to kaons. All these simplifications 
yield a somewhat conservative result for the conventional AN spectra and must be avoided in the future study; the 
proper inclusion of the kaon regeneration is the most essential point. 

(v) At the stage of nuclear cascade (but of course not at the lepton production stage) the decay of charged pious 
is neglected. This approximation greatly simplifies the description of the pion regeneration/overcharging and nucleon 
production in pion-nucleus collisions and it is valid at the pion energies above 1-2 TeV for directions close to vertical. 
However it becomes too crude for near-horizontal directions at Ei, < 7 — 8 TeV. To extend our results up to 1 TeV, 
we included the appropriate corrections for pion decay using a numerical procedure. 



V. AN FLUX (NUMERICAL RESULTS) 

Let us discuss the numerical results presented in figures |^-^ 

Figures!^ and ^ show the individual contributions from n^2, K^2i ^"31 ^"^'^ f^ts decays into the Ve+T^e and I'^+Fp 
fiuxes for vertical and horizontal directions, as well as the overall tt, if-neutrino fluxes. As an example of possible FN 
contribution, we also show the results obtained with the three alternative models for charm hadroproduction: the 
recombination quark-parton model (RQPM), quark-gluon string model (QGSM) and the model based on perturbative 
quantum chromodynamics (pQCD). The basic assumptions of the first two models were described by Bugaev et al. [7] 
(see also and references therein). The third, "state-of-the-art" model was proposed recently by Thunman et al. \i] 
to simulate charm hadroproduction through pQCD processes. To leading order in the coupling constant, as, these 
are the gluon-gluon fusion and the quark-antiquark annihilation. The next-to-leading order contributions are taken 
into account by doubling the cross sections. To simulate the primary and cascade interactions, the authors use the 
well-accepted Monte Carlo code PYTHYA. 

It is our opinion that the RQPM and QGSM give the safe upper and lower limits for the prompt muon and 
neutrino fiuxes. These limits are not inconsistent with the current deep underground measurements of the muon 
intensity. However, considering rather strong discrepancy between the data of the ground-based and underground 
muon experiments 0], the comparatively low prompt muon contribution predicted by the pQCD model cannot be 
excluded. Similar, very low prompt muon contribution has been evaluated recently by Battistoni et al. Q using the 
DPMJET-II code based on the two-component dual parton model and interfaced to the shower code HEMAS. The 
calculation of Battistoni et al. |^ does not yield the absolute prompt muon flux but, from the estimated prompt- 
to-conventional muon ratio, one can see a leastwise qualitative agreement with the result of the pQCD model by 
Thunman et al. Q. In particular, both models predict that the prompt muon contribution overcomes the vertical 
TT, K-imion fiux in the region of a thousand TeV and therefore is undistinguished in the present-day muon experiments. 
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As may be seen from fig. |[ the K^s decays give the main contribution into the conventional + Ve flux above 
1 TeV independent of zenith angle. The if^ and K'^-^ contributions are practically equal at i9 = 0° and close in 
magnitude at i? = 90°, despite the 8-fold difference between the and if^g decay branching ratios. The first reason 
is that the life time of K\ is about 4.2 times that of . The second one lies in the different inclusive cross sections 
for + K~ and if" + K'^ production (see Vail et al. [|l6|). For the conventional v^^ flux, the main contribution 
comes from decays (fig. ij). At = 90° it stands out above the PN contribution predicted by RQPM, up to about 
200 TeV. Second in importance is the 7r^2 decay contribution. The cross sections for pion production are almost order 
of magnitude larger than for kaon production, but pion decays become rare above f TeV owing to the large life time 
and Lorentz factor compared to kaon ones. The iC^a decays give comparatively small contribution which however 
cannot be neglected. Muon decay contribution is negligible at = 0° and very small at ?? = 90°. 

In figures ^ (a-d) we present the conventional i^e + and //^ + 77^ fluxes at = 0° and 90° as calculated by 
Volkova H], Mitsui et al. ||^, Butkevich et al. Q], Lipari |^ and Agrawal et al. All these fluxes are normalized 
to our one. In the multi-TeV energy range and above, our results fall within the lowest and highest predictions. At 
Ey < 2 — Z TeV and < 10 TeV for, respectively, d — 0° and — 90°, our calculations give somewhat excessive 
i^ij. + 1^11 fluxes. In part, this is due to our simplified consideration of the pion decay effect (see Section ^ ). The 
corresponding discrepancy is almost negligible for electron neutrinos, since pions give no direct contribution to the 
i^e + i^e production. 

In fig. 1^ (a-d) we show our result for the conventional neutrino to antineutrino ratios at i? = 0° and 90° (vs. 
neutrino energy) together with the results of Butkevich et al. Q and Lipari 1^. The Ve/T^e and v^/v^ ratios are very 
sensitive to the input parameters like the n/p ratio, which is governed by the primary chemical composition, and the 
tt^/tt", /K~ , K^/K^ ratios, which are determined by the cross sections for the meson production, regeneration 
and overcharging. Whilst the parameters adopted by Butkevich et al. |Q, Lipari |^ and in the present calculation 
are all within the limits dictated by the available cosmic-ray and accelerator data, the sets of the parameters differ 
considerably from each other. Consequently, it is no great surprise that the predictions of different models vary over 
a wide range. 

On the other hand, owing to the well-known difference between the cross sections for i^iN and ViN charged and 
neutral current induced interactions, the Ve/i^e and Vfi/V^ ratios are very important inputs for the correct evaluating 
the neutrino-induced throughgoing muon flux and contained event rate in the neutrino detector ^ As is seen from 
fig. 1^, the present-day uncertainty in the v^/V^ ratio is not satisfactory. The situation with the Ve/Vf. ratio proves to 
be somewhat better. 

Figures ^ (a,b) show the neutrino flavor ratio, i?, vs energy at t9 = 0° and 90°. The dashed curves represent the 
ratio for the conventional AN flux which is a monotonically increasing function of energy varying from about 27.6 to 
34.3 at I? = 0° and from about 13.3 to 33.6 at = 90° within the interval 1 100 TeV. The sohd curves are for 
the R evaluated with taking into account the PN contribution according to three abovementioned models for charm 
production. The PN contribution results in a decrease of the AN flavor ratio, because the semileptonic decay modes 
of charmed hadrons with (V^) and Ve (T^e) in the final state have almost the same branching ratios. 

As one can see, the R is very sensitive to the charm production model even at energies where the PN contribution 
remains small in comparison with the conventional one. This effect provides an interesting potential possibility for 
the experimental discrimination of the PN production by measuring the ratio of the "muonfuU" (t'^ and induced) 
to "muonless" (j/g and Ve induced) contained event rates in a neutrino telescope for different energy thresholds and 
directions. To a certain extent, such an experiment is easier than the measurement of the absolute neutrino event 
rate. 

Let us briefly sketch now the Kg^ form factor effect. Our calculations show that the effect under discussion is almost 
independent of zenith angle but quite different for muon and electron neutrinos. At E^, > 1 TeV, the effect is identical 
for the K^2, ^'Ud contributions as well as for the overall + Ve flux (since the Kes, decays are the main source of 
Ve and Fe); its magnitude is just higher than 3%. The magnitude of the effect is different for K^^^ (about 6%) and 
ii'^g (about 4%). But, considering that the iiT^g contribution is by itself small, there is no any change in the overall 

+ flux. As a consequence, the inclusion of the g^-dependent Kg'^ form factors decreases the {v^ +F^)/(j^e +Fe) 
ratio for the conventional neutrinos by about 3 to 4%, depending on the zenith angle and energy. As one might 
expect, the effect is small but not completely negligible. Figure |g illustrates the Kg^ form factor effect for the AN 
flavor ratio in the presence of a PN contribution. Clearly, it is absent when prompt neutrinos dominate. Thus, we use 
the lowest PM contribution as predicted in the pQCD model by Thunman et al. By the evident reasons, there is 
no effect of any value in the Ve/T^e and v^jv^ ratios. 



^At PeV energies, the Ve/i'e ratio becomes a crucial parameter on account of the W resonance in Uee annihilation. 
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VI. CONCLUSIONS 



The main result of this work consists in the expUcit formulas for the semileptonic decay rates (differential and 

total) of and Kl mesons, which take into account the g^-depcndent Kf^ form factors. The obtained rates differ 
essentially from those calculated with constant form factors and the total rates are in good agreement with experiment. 

With a rather detailed model for nuclear cascade in the atmosphere and with the improved differential K^^ rates, 
we have calculated the i/e, Fg, i^^ and 77^ spectra at energies 1 to 100 TeV. The calculated spectra are within the 
limits resulting from the uncertainties in the current data on the primary cosmic ray flux and cross sections for 
hadron-nucleus interactions at high energies. 

The outcome of the inclusion of g^-dependent K^^ form factors into the AN flux calculation is as follows. The 
electron neutrino flux from K% and Kf^ decays increases by about 3-3.5%. The K^^ and -FT^g decay contributions 
into the muon neutrino flux reduces by about 6 and 4 %, respectively, whereas the change in the overall ly^ -\-V^ flux is 
completely negligible. If the PN contribution is as slight as predicted by the pQCD model, the change in the neutrino 
flavor ratio, i?, comprises 3 to 4%. This small systematic effect may only slightly be enhanced by taking account of 
kaon regeneration or through the variation of the kaon production cross sections within the experimental boundaries. 
And vise versa, it may be removed, wholly or in part, beyond the multi-TeV energy range if the PN contribution is 
as large as it follows from the RQPM or QGSM predictions. 

The sensitivity of the R to the PN contribution provides a way for an experimental discrimination of this contribution 
by measuring the relationship between the "muonfuU" and "muonless" contained event rates in a neutrino telescope. 
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